Spin-1 Kitaev model in one dimension 
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We study a one-dimensional version of the Kitaev model on a ring of size N, in which there is a 
spin S > 1/2 on each site and the Hamiltonian is JJ] n S^S^ +1 . The cases where S is integer and 
half-odd-integer are qualitatively different. We show that there is a Z2 valued conserved quantity 
W n for each bond (n, n + 1) of the system. For integer S, the Hilbert space can be decomposed 
into 2^ sectors, of unequal sizes. The number of states in most of the sectors grows as d , where d 
depends on the sector. The largest sector contains the ground state, and for this sector, for S = 1, 
d = (V5 + l)/2. We carry out exact diagonalization for small systems. The extrapolation of our 
results to large N indicates that the energy gap remains finite in this limit. In the ground state 
sector, the system can be mapped to a spin- 1/2 model. We develop variational wave functions to 
study the lowest energy states in the ground state and other sectors. The first excited state of 
the system is the lowest energy state of a different sector and we estimate its excitation energy. 
We consider a more general Hamiltonian, adding a term A J~] W n , and show that this has gapless 
excitations in the range Ai < A < A§- We use the variational wave functions to study how the 
ground state energy and the defect density vary near the two critical points Ai and A2. 



PACS numbers: 75.10.Jm 



I. INTRODUCTION 



In recent years, there have been many studies of quan- 
tum spin systems which are characterized by a high de- 
gree of frustration and topological order. The word 'frus- 
tration' here refers to systems with competing interac- 
tions having a large number of states with energy near 
the minimum energy. Topological order implies the ex- 
istence of invariants which, for topological reasons, are 
robust against a large class of perturbations. Such sys- 
tems are often associated with a novel structure of the 
ground state and low-lying excitations, and are interest- 
ing from the point of view of possible applications in 
quantum computation A particularly interesting 

model in this context is the two-dimensional frustrated 
spin-1/2 model introduced by Kitaev |3j. This model 
has several fascinating properties which have been stud- 
ied in great detail 12} . For instance, the model and its 
variants constitute the only known class of spin models 
in two dimensions or more dimensions that is fully in- 
tegrable, being reducible to a system of non-interacting 
Majorana fermions. A similar model, called the compass 
model, although not exactly solvable, was introduced by 
Kugel and Khomskii many years ago 13] to understand 
the magnetic properties of transition metal oxides which 
have orbital degeneracies. Recently physical realizations 
of the spin-1/2 Kitaev model have been proposed in op- 
tical lattice systems [l4j and in quantum circuits [in- 
variants of the model have also been studied in two di- 
mensions [l6|-[24j , three dimensions pol [26j | and also on 
quasi-one-dimensional lattices |27H29l |. Finally, the spin- 
s' Kitaev model has been studied in the large S limit 
using spin wave theory [30(, and the classical version of 
the Kitaev model has been studied at finite temperatures 
using analytical and Monte Carlo techniques |31|. Their 



results indicate that while the phenomenon of order-by- 
disorder [32H3a ] may occur in the quantum mechanical 
Kitaev model, it does not in the corresponding classical 
model. 

For the Kitaev model with spin S > 1/2, there is a 
Z2 invariant associated with each plaquette for arbitrary 
spin-S, which reduces to the conserved Z2 gauge flux 
for the spin-1/2 case [3(| ■ However, the model does not 
seem to be fully integrable. While some differences in 
the structure of the invariants between the models with 
half-odd-integer and integer spins have been pointed out 
[30j , the issue of whether there are systematic differences 
in the nature of the low-energy spectrum is also of in- 
terest. In the present paper, we approach this prob- 
lem by examining the spin-1 Kitaev model. The two- 
dimensional model appears difficult to analyze, but even 
the one-dimensional version of it has a lot of interesting 
structure, as we proceed to show. 

The plan of this paper is as follows. In Sec. UH we 
consider the spin-S Kitaev chain. In Sec. Ill AI we show 
that this model has local, mutually commuting conserved 
quantities W n , for integer S. The eigenvalues of W n are 
±1. For open boundary conditions, there are some ad- 
ditional conserved quantities at the ends of the system. 
The existence of these conserved quantities implies that 
the Hilbert space of a JV-site system can be decomposed 
into a sum of 2 N disjoint subspaces. The dimensions of 
these subspaces are not equal. In Sec. Ill B I we develop 
a formalism to compute the dimension of these sectors. 
For large N, the dimension varies as d N in most sectors, 
with the constant d depending on the sector. The sec- 
tors show complicated spatial structures, arising from the 
spatial structure of {W n }- We show this in Sec. Ill CI by 
computing the non-trivial spatial dependence of expec- 
tation values of spin operators in some sectors, averaged 
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over all states in the sector. We then consider the spin-1 
model in Sec. IIIII In Sec. IIII A[ we consider the ground 
state and lowest excited state of the system. Exact diag- 
onalizations of small systems show that the ground state 
lies in a sector in which W n = +1 for all n. In this sec- 
tor, there is a gap between the ground state and the first 
excited state. The lowest excited state of the system is 
the ground state of a different sector; and the energy- 
gap seems to approach a non-zero value in the limit of 
the system size going to infinity. In Sec. IIV[ we con- 
sider the sector containing the ground state, and show 
that the Hamiltonian is equivalent to the Hamiltonian of 
a deposition-evaporation process of a nearest-neighbor- 
exclusion lattice gas model, which can be written as of a 
spin-I/2 system with local interactions with a range ex- 
tending to at most next-nearest neighbors. The Hamil- 
tonian seems to be difficult to diagonalize exactly, we 
present a variational study of the ground state in Sec. [V] 
The variational estimate of the ground state energy is 
found to agree well with the results obtained numerically 
for small systems. We also analyze the first excited state 
of the Hamiltonian. In Sec. IVI1 we consider a more gen- 
eral Hamiltonian, obtained by adding a term AJ^_ W n , 
and discuss its ground states as a function of A. We show 
that the ground state of this new Hamiltonian is gapless 
for a range of couplings A^j; < A < A§, and gapped oth- 
erwise. We argue that for A just above AJ, in the sector 
containing the ground state, the density of negative Ws 
is of order I i g(A~A c ) I' ^ or ^ J us ^ below Aj, the density 
of positive Ws goes to zero as (A| — A) 1 / 2 . In the final 
section, we summarize our conclusions, and discuss the 
relationship of this model with the Fibonacci chain. 



II. ONE-DIMENSIONAL KITAEV MODEL 

In this section, we will discuss a one-dimensional spin- 
s' model which is obtained by considering a single row of 
the Kitaev model in two dimensions. 

Let us begin with the Kitaev model on the honeycomb 
lattice. This is governed by the Hamiltonian 
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where (ij) a denote the nearest-neighbor bonds in the a th 
direction. If we set J z = 0, we get a set of decoupled 
chains. We call this the Kitaev chain, and this is the 
topic of this paper. The Hamiltonian is by 
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In general, the couplings J m could be all different from 
each other. If some of the couplings are negative, we can 
change the signs of those couplings by performing the 
unitary transformation 
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FIG. 1: Picture of the Kitaev chain showing one of the con- 
served quantities Wn- 



on appropriate sites. We consider the simpler case, where 
all couplings have the same value, J m = J. Without any 
loss of generality, we set J = 1. Finally, the Hamiltonian 
can be unitarily transformed to a more convenient form 
by the following transformation on the even sites, 



QX . QV 
°2n ~ f °2m 



(4) 



The Hamiltonian in Eq. ([2]) then takes the translation 
invariant form 



(5) 



A. Invariants 

The Hamiltonian in Eq. (JS|) has the following local 
symmetries for all S. Let us introduce the operators on 
sites 



— e I 

and operators on bonds 

w — y.v y x 

as shown in Fig. 1. We then find that 
[W n ,H} = 0. 



(6) 
(7) 
(8) 



The eigenvalues of S° are ±1 for integer S and ±i for 
half-odd-integer S. Thus for any value of the spin S, the 
eigenvalues of W n are ±1. 

However, there is a qualitative difference between inte- 
ger and half-odd-integer values of S. For integer values of 
S, all the matrices S° matrices commute with each other, 
whereas for half-odd- integer values, commutes with 
Yj b m for n 5^ m but anticommutes with for a ^ b. Con- 
sequently, for integer S, all the invariants W n commute, 
but for half-odd-integer S, W n anti-commutes with its 
neighboring invariants, W n ±i, and commutes with W m , 
m =/= n, n± 1. We will now show that this implies that all 
the eigenstates of the chain with half-odd-integer S are 
2 N / 2 fold degenerate. 

The invariants for half-odd-integer S can be combined 
in the following way to form a set of mutually commuting 
angular momentum operators, one per every two bonds, 



< = W 2n , ill = W 2n _i Y[ Wj 



2m- 1, 



( 9 ) 
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It can be verified that 



= 2*5. 



06c c 



(10) 

(ii) 



The /i^ commute with the Hamiltonian as they are made 
by multiplying conserved operators. Hence Eq. (fTU|) 
shows that the Hamiltonian has a (SU(2)) N ^ 2 symme- 
try, where N is the number of sites. Eq. (fTTj) shows that 
each of these SU(2) factors are realized in the spin-1/2 
representation. Thus each eigenstate is 2 JV / 2 -fold degen- 
erate. There is no reason for such a degeneracy for integer 
S and indeed, as we will see later, the ground state for 
S = 1 is non-degenerate. 

We note that the spin-S* Kitaev model in two dimen- 
sions also has a Z 2 valued invariant associated with every 
hexagon of the honeycomb lattice [3(| ■ When they are 
restricted to a single chain, the invariants take the form 



v — v9y z v 



+ l^n+2 



(12) 



which involves three neighboring sites. The invariants 
given in Eq. (J7J are simpler because they only involve 
two sites. For any spin S, we find that E^E^E* = I and 
(E^) 2 = (— 1) 2S ; hence the invariants in Eqs. J7]) and 
(fT2")l are related to each other as 



V n = (-1) 2S W n W n+1 



(13) 



Open chains have some extra symmetries at the edges. 
If the site labels of the open chain are 1, • • • ,N, then Sf 
and S N also commute with the Hamiltonian. Thus at 
the first and last sites, we have a U(l) symmetry group 
generated by these operators. Note that a Z 2 subgroup 
of this group, consisting of the operators Ef and E^, 
also commutes with all the invariants. If we combine the 
operators Sf and S N with the W n invariants on the first 
and last bonds, we have a larger symmetry group made 
of Wi, Sf and their products at the first bond, and the 
group made of Wn, S n at the last bond. As we will 
show in Sec. IIII1 for the 5 = 1 case the group formed is 
517(2) x U{1) at each end. 



B. Counting of states for integer S 

We will now develop a formalism to count the number 
of states in a given sector for integer S. In this case, the 
E matrices commute and hence can be simultaneously 
diagonalized. If \S, m) denote the eigenstates of S z , then 



Y?\S,m) = (~l) s \S~m), 
Z y \S,m) = (-l) s+m \S,-m) 
Z z \S,m) = (-l) m \S,m). 



(14) 
(15) 
(16) 



We can construct the eigenstates of the E matrices in the 
m =/= subspace, 



\S,m±) = {\S,m) ± \S,-m)) 



(17) 



where m = 1, • • ■ ,S. The eigenvalues of the matrices are, 

Tf\S,m±) = ±(-l) s \S,m±), (18) 
H v \S,m±) = ±{-l) s+m \S,m±), (19) 
E z |S*,m±) = {-l) m \S,m±). (20) 

The states of a chain can be classified by the eigenval- 
ues of Eif and EH as 



(VnXn) ■ ■ ■ (y2X2){yiXi), 



(21) 



where x n ,y n — ±1 are the eigenvalues of E^ and T, n 
respectively. The invariants are then W n = x n+ iy n . 

We now calculate the number of states in a given sector 
W using a standard transfer matrix technique. Consider 
the allowed states of r sites, when the values of r — 1 
constants Wj , with j — 1 to r— 1 have been specified. We 
denote this set of values by W. Let Z r (y\W) denote the 
number of allowed states of this set of sites with Ejf = y 
where y takes values ±1. We now add a site r + 1 to the 
chain, and also specify W r . Let the new set of {W} be 
denoted by W. 

Consider first the case W r — +1. Clearly, we can have 
two possibilities: E^ +1 = EJf = +1, or E 



r+l 



EH 



1. 



Let v{p,p') denote the number of states of a single site 



with £« = p, and E* 
recursion equation 



p' . Then, we clearly have the 



Zr+i{y\W) = v{y, +l)Z r (+l, W) + v{y, -l)Z r (-l, W). 
This equation can be written as a matrix equation 



z r+1 (+i\w) 
z r+1 {-i\w) 



T 4 



Z r {+l\W) 

Z r (-l\W) 



where T + is a 2 x 2 matrix given by 

T + = 



K-1,+1) K-1,-1) 



It then follows from Eqs. (|18ti20|) that 



T+ = - 
+ 2 

1 

2 



S-l S- 
S + l s- 

s + 
s 



s 
s 



for S odd, 
for S even. 



(23) 

(24) 

(25) 
(26) 



Similarly, when W r = 
equation is 

z r+1 (+i\w) 

Z r+1 (-l\W) 



T, the corresponding recursion 

(27) 



Z r (+l\W) 
Z r (-l\W) 



where the matrix T_ is given by 
T-=T>t x , with 



1 

1 



(28) 



It is then clear that for a given set of invariants W, the 
number of states can be written in terms of a product of 
the matrices T+ and T_. 
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For example, for an open chain of N sites, and W = 
{W N _ 1 ,...W 3 ,W 2 ,W 1 } = { + 1, ... + 1, -1, -1}, we have 



Z N (-1\ 



T+...T+T-T- 



Zi(+l\<l>) 
Zi(-l|0) 



(29) 



where <fi denotes the null string, and Zi(y\(f)) denotes the 
number of states of the spin at site 1 with Yi\ = y. Thus 
Zi(+l|0) = S+l, Zi(-l\cf>) = S, when S is an even 
integer, and Zi(+l|^) = S, Zi(-l|0) = S + 1 when S is 
an odd integer. The total number of states in this sector 
is then given by 



r(w) - z N (+i\w) + z N (-i\w). 



(30) 



For a closed chain, there is an additional invariant 
Wn — VnX\ and the number of states in the sector be- 
comes 



N 



T(W) = Tr m T 



w„ 



(31) 



where fw n = T± for W n = ±1 and I]_=i is an ordered 
product of T± matrices, from site 1 to N with the index 
increasing from right to left. 

We now calculate the dimensions of some sectors for 
a closed chain of length N. It is easy to get an explicit 
answer for the two extreme limits when W n — ±1 for all 
n. In these cases, the number of states, T^, is 



r± = (<£)" + 



(32) 



where df(S) and d% (S) are the larger and smaller eigen- 
values of T± respectively. The eigenvalues can be com- 
puted to give, 



U l(2) 



l l(2) 



- (S ± VS 2 + 2S + 2) for S odd, (33) 

- (s + 1 ± VS 2 + l) for S* even, (34) 

- (S + 1 ± v 7 ^ 2 - l) for 5 odd, (35) 

- (s ± VS 2 + 2S*) for S even. (36) 



2 
1 

2 
1 

2 
1 

2 

For S = 1, df is equal to the golden ratio, 7 = (1 + 
VS)/2, and 4 = -I/7. As iV — > 00, the dimension of 
the Hilbert space in the sector with all W n = 1 grows as 
j N . On the other hand, = d^ = 1. The dimension of 
the sector with all W n = — 1 is therefore equal to 2. 

With the exception of S = 1, the larger of the two 
eigenvalues df 1 is always greater than 1, and in the N — > 
00 limit, we have 



r±(s) = (df(s)) 



N 



(37) 



d l {S) is referred to as the quantum dimension of the 
sector. As can be seen it is, in general, fractional for any 
S. In the limit S — > 00, the quantum dimension tends 
to S + 1/2 for both the sectors. It is interesting to note 
that it is a half-odd-intcgcr in this limit. 



Expectation values of the E operators in 
different sectors 



In this section we find the expectation values of the 
operators in various sectors. We will assume peri- 
odic boundary conditions. Our calculation will average 
over all the states of a given sector considered with equal 
weight; this can be considered as a calculation in the limit 
that the temperature T — > 00, so that it does not depend 
on the Hamiltonian. 

We evaluate the expectation values of by inserting 
projection operators at site n in the product of transfer 
matrices in Eq. pip . This yields the following expression 
for the expectation value of the E° operator in a general 
sector with a VK-configuration W 



w 



Tr 



n t w, ] ^„ ( n j ^ 

, j=n+l 



i=l 



/r({w}) 



where 

T x Wn = W n T Wn T* 
T v Wn = r z T Wn 



T 



W n T Z T Wn T Z 



(38) 



and t z and t x are the well-known Pauli matrices. 

We now compute the expectation values of in two 
sectors: the sector Wo with all W n — +1, and the sector 
Wi in which one of the W n = — 1 and all the other W n 
■ I (without loss of generality we pick Wn = — 1)- The 
expressions for (E°)yy 1 = (E°)o,i can be evaluated in 
terms of the eigenvectors and eigenvalues of T + . The 
T + matrix is a linear combination of the Pauli matrices, 
t z and t x . Its eigenvectors are spinors polarized parallel 
and anti-parallel to a direction in the z x plane, forming 
an angle 9$ with the z— axis, where #s is defined by, 



cos t>s = — 



1 



for S even, 



for S odd, 



sm t>s 



Vl + S 2 
5 + 1 



s 



for S odd, 



VTTs 2 



for S even. 



(39) 



(40) 



For the sector with all W n — +1 it is easy to see that 
(E*) = (T,y) . For large N, we obtain 



(E^ } ) 



COS 6*5, 



(41) 
(42) 



(E*) = cos^ 0s + %bw*6 s . 
In the sector where Wn = — 1 and the rest are equal 



5 



A 
(0 c 



0.8 



V 

0.4 
0.2 




+ + + + + + + - + - 

'1 3 5 7 9 11 13 15 17 

n 



FIG. 2: Plot of (S* 2 ) (dotted line), (S£ 2 ) (dashed line) and 
(S^ 2 ) (full line) as a function of n for S 1 = 1, on a ring with 

16 sites in the sector + + + H 1 h H 1 — , with 

periodic boundary conditions (site 17 = site 1). 



to +1, we get, for large N, 



(K)i 



1- 



dt 



2 cos 9 s 



d{ cosg +d% sin 2 S \ d 



d+ 



Note that in the limit n -> ±oo, (£") in Eqs. (|43H45|) 
approach the values given in Eqs. (|4HI42[) exponentially 
quickly. 

While in general S° are complicated multi-spin opera- 
tors, for S = 1 we have ££ = 1- 2(S%) 2 . Thus, for S = 1 
we are essentially computing the expectation values of 
(S 1 ^) 2 . To see what the spin textures are like in a typical 
sector, we have plotted in Fig. [3] the expectation values 
of S% 2 , S% 2 and S* 2 for S = 1, as a function of the spatial 
coordinate n for a ring of size 16, in the sector where the 

sequence of W's is H — I — I — I 1 hH 1 — . 

This sequence was chosen as it is a dc Bruijn sequence 
J>6j of length 16, in which each of the 16 possible bi- 
nary sequences of length 4 occur exactly once, taking the 
periodic boundary conditions into account. 



III. 



1 MODEL 



We will now focus on the Kitaev chain with spin-l's at 
each site. We will work with the natural spin-1 represen- 
tation in which 



(46) 



In this representation, the matrices E° are diagonal and 
are given by 




(47) 



We note that these matrices satisfy E x E y E z = /. Wc 
denote the basis vectors by |x), \y) and \z) defined as 



I*) 



(48) 



We then see that the 9 possible states at sites (n, n+X) 
are given by: 

\xy), \xz), \yx), \zy) and \zz) with W n = 1, (49) 
and 

\xx), \yy), \yz) and \zx) with W n = — 1. (50) 

From Eq. (|4T>)) we have, 

S*\x) = 0, Sy\x)=i\z), S z \x) = -i\y), 

S x \y) = -i\z), S»\y)=0, S*\y)=i\x), 
S x \z) = i\y), S*\z) = -i\x), S*\z)=0. (51) 

Eqs. (gll) and $5T$ imply that (S a ) 2 = (1 - E a )/2. 

For the 5 states in Eq. (|49j) satisfying W n — 1, we 
have the following actions of the relevant term in the 
Hamiltonian, 

S?S$\xy) = 0, 
SfSl\xz) = 0, 
SfS%\zy) = 0, 
SfS$\zz) = \yx), 

SlS y 2 \yx) = \zz). (52) 



For the 4 states in Eq. ([50|) satisfying W„ = —1, the 
actions of the relevant term in the Hamiltonian are given 
by 

SfS$\xx) = 0, 
S?S$\yy) = 0, 
SfSi\yz) = -\zx), 

= -\yz). (53) 



SfS y 2 \zx) 



As mentioned earlier, for an open chain with site num- 
bers going from 1 to N, we hnd that Sf and 5^ commute 
with H . We dehne the operators, 



■SfWiSf. 



EE Sf. 



(54) 



r°EE-(l-(Sf) 2 ). (55) 
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It can be verified that these operators obey a SU(2) x 
U(l) algebra. Exactly the same construction on the last 
bond, with Sf —> S N and W\ — > Wn, yields the same 
algebra on that bond. 



A. Numerical studies 

We have carried out exact diagonalization studies of 
small systems with periodic boundary conditions in or- 
der to find the energies of the ground state and the lowest 
excited state of the spin-1 Kitaev chain. We find that the 
ground state lies in the sector with all W n = 1 and has 
zero momentum (momentum is a good quantum number 
in this sector since the values of the W n 's are translation 
invariant). The ground state energy per site as a func- 
tion of the system size N is presented in Table I. We see 
that Eq/N shows odd-even oscillations as a function of 
N but seems to converge quite fast. The fast convergence 
indicates that the ground state must have a fairly short 
correlation length. The TV-dependence of En = Eq/N 
can be fitted to the form 



E N =E oc + B(-a) 



N 



(56) 



A simple plot of log \ En — Eoo | versus N (Fig. [3]), gives a 
good straight line for E^ = —0.60356058, which we take 
to be the best estimate of E^, . The corresponding values 
of B and a are 0.07 and 0.51. The estimated errors of 
extrapolation are about 1 in the last significant digit. 



N 


Eo/N 


N 


Eo/N 


2 


-0.707106770 


11 


-0.603525102 


3 


-0.577350259 


12 


-0.603578389 


4 


-0.612372458 


13 


-0.603551567 


5 


-0.600000024 


14 


-0.603565216 


6 


-0.605046094 


15 


-0.603558183 


7 


-0.602888465 


16 


-0.603561819 


8 


-0.603869855 


17 


-0.603559971 


9 


-0.603412688 


18 


-0.603560924 


10 


-0.603632331 


19 


-0.603560388 



TABLE I: Ground state energy per site versus N. 

In the sector with all W n = 1, the first excited state 
has momentum equal to ir if N is even. We find that 
the gap separating it from the ground state is given by 
1.0353 for N = 4 and 0.9845 for N = 6. These values 
also seem to be converging rapidly, and the large value is 
consistent with a short correlation length. However, this 
is not the lowest excited state of the system. Rather, 
we find that the state nearest in energy to the ground 
state is the ground state of the sector with exactly one 
W n = — 1 and all the other W n — 1. (We cannot use 
momentum to classify the states in this sector since it is 
not translation invariant). The energy gap AE between 
the lowest energy state in this sector and the ground state 
of the sector with all W n = 1 is shown in Table II. We 




FIG. 3: Graph of \E /N - E^ with N, where E^ = 
—0.60356058, showing an exponential convergence to the 
value of the ground state energy with N. 



see that these also oscillate between even and odd values 
of N but seem to converge quite fast to a small but non- 
zero value. This is evidence that the spin-1 Kitaev chain 
has a finite gap in the thermodynamic limit N — > oo. 



N 


AE 


3 


0.1141 


4 


0.2025 


5 


0.1671 


6 


0.1802 



TABLE II: Energy gap between ground state and first excited 
state versus N. 



IV. MAPPING THE SPIN-1 CHAIN TO A 
SPIN-1/2 CHAIN 

For a given value of the state of the spin at site n, 
and a given value of W n , there are at most two choices 
for the spin state at site n + 1. Hence it is clear that 
the Hilbert space of a given sector can be mapped into 
the Hilbert space of a spin- 1/2 chain, with some states 
excluded which correspond to infinite energy. However, 
in general, the corresponding Hamiltonian would have a 
rather complicated form, with long-ranged interactions. 
The mapping is easy to construct explicitly in the sector 
with all W n = +1, and the corresponding Hamiltonian 
has only local interactions. This is what we now proceed 
to show. 

Consider the state zzzz ■ ■ ■ that belongs to the sector 
with all W n = +1. The only allowed process in this sector 
is zz ^ yx [Eqs. (|52p]. We may think of this process as 
a quantum dimer deposition-evaporation model. The z- 
spins are treated as empty sites; two empty sites can be 
changed to being occupied by a dimer yx by a 'deposition' 
process, and conversely, yx can 'evaporate' and become 
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zz again. The dimers have a hard-core constraint, and 
a site cannot be shared by two dimers. The dimers are 
oriented: the 'head' x being to the right of the 'tail' y. 

This dimer deposition-evaporation model can also 
be described as a deposition-evaporation of a nearest- 
neighbor exclusion lattice gas. We just think of the heads 
as particles, and do not distinguish between the tails and 
empty sites, except for ensuring that we deposit a par- 
ticle at a site only if it is empty and both its nearest 
neighbors are also empty. Then this model is described 
by the Hamiltonian 



Ha 



7 E 



(57) 



We note that this model is different from the dimer 
deposition-evaporation models studied earlier (37| . in 
that the two ends of the dimer are distinct, and there is 
no reconstitution. Also, this Hamiltonian does not have 
an interpretation as the evolution operator of a classical 
Markov process, as there are no diagonal terms corre- 
sponding to probability conservation. 

We have introduced a minus sign in the Hamiltonian 
for later convenience. This does not change the eigen- 
value spectrum as the eigenvalues of Hd occur in pairs 
±ei. 



V. VARIATIONAL STUDY OF SECTOR WITH 

ALL W n = 1 

We will now use a variational approach to study the 
ground state of the Hamiltonian Hd with periodic bound- 
ary conditions. We use the z-basis, and denote the f state 
at the site i by an occupied site (rij = 1), and the J, state 
by an empty state (n, = 0). Since two adjacent sites can- 
not be simultaneously occupied, the state space is that 
of hard-core particles with nearest-neighbor exclusion on 
a line. A configuration C is specified by an iV-bit binary 
string 0010010101 • • • , which gives the values of all the N 
occupation numbers rij. We note that in the basis where 
all the rii are diagonal, the Hamiltonian Hd has all ma- 
trix elements non-positive. This implies that the (real) 
eigenvector corresponding to the lowest energy will have 
all components of the same sign in this basis. 

For the ground state of Hd, we consider a variational 
wave function of the form 



£v^ob(C) |C) 



(58) 



where Prob(C) is chosen as the probability of the lattice 
gas configuration C in some classical equilibrium ensem- 
ble corresponding to a suitably chosen lattice gas Hamil- 
tonian. Clearly, this trial vector is normalized, with 



(59) 



With this choice, Prob(C) is also the probability of the 
configuration C in the quantum mechanical variational 
state lib). 



The simplest choice of the lattice-gas Hamiltonian is 
that of a classical lattice gas with nearest-neighbor ex- 
clusion, and a chemical potential /i, with a Hamiltonian 
given by 



H r l = 



-co 2J nin i+ i 



(60) 



where we use the convention that • oo = 0; hence the 
first term in Eq. (ISTfl) allows states with n^ni+i = 
but disallows states with mrii+i = 1. Let us denote 
z — exp(/3/x). It is straightforward to determine various 
correlation functions in the thermal equilibrium state cor- 
responding to H c i . The probability of a configuration C 
is given by 



Prob(C) = eM~PH cl {C)]/n N {z) 



(61) 



where f2jv(z) is the grand partition function for a ring of 
N sites. 

The grand partition function Qn(z) can be determined 
using the standard transfer matrix technique. We find 
the largest eigenvalue of the 2x2 matrix T2 given by 



1 1 

z 



(62) 



We now calculate (ip\Hd\ip). The matrix element of the 
i-th term is clearly zero, unless n,_i = nj+i = 0. Then 
the only non-zero matrix element is 

(H d ) /N = -%Jz Prob(000) = ^Prob(OlO). (63) 



Here Prob(OOO) denotes the probability that randomly 
selected three consecutive sites in the ring will be empty 
in the classical ensemble, and similar definition for 
Prob(010). This is easily calculated for the Hamiltonian 
H c i in the limit of large N. We get 



Prob(010) = Prob(l) = p. 
The largest eigenvalue A of T2 is given by 



A = (1 + VT+4i)/2, 



(64) 



(65) 



and p is the density per site given by p = zd\og(A)/dz. 
Extremizing (Hd) with respect to z, we find that the 
minimizing value occurs for z = 0.405, yielding (Hd) = 
—0.60057. This gives us the variational bound variational 
bound on the ground state energy per site Eq 



E a < -0.60057. 



(66) 



This energy is somewhat higher than the energy obtained 
in the previous section (see Fig. [3]), indicating that the 
correlations in the classical Hamiltonian H c \ do not ex- 
actly reproduce the correlations in the quantum ground 
state of Hd- 
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We can make a better variational calculation by con- 
sidering a classical lattice gas with an additional next- 
nearest-neighbor interaction. The Hamiltonian of this 
lattice gas is 

H' cl = +oo^n;7i m - K^Tiini+2 - fJL^^m. (67) 



eigenvector in the sector with all W n — +1. In both 
these regions, the gap in the excitation spectrum is of 
order |A|. As we vary A from — oo to +oo, initially the 
gap decreases and becomes zero at some value Af . We 
then expect a gap to open up again when A is greater 
than a second critical point \? 2 > XI ■ 



Let us denote z = exp(/3/i), and u = exp(/3K). In this 
case, the transfer matrix is a 3 x 3 matrix given by 



1 1 

z zu 
1 



(68) 



The probability of the configuration C in the equilibrium 
ensemble is given by 



Prob(C) = eM-0Kl(C)}/n N (z,u), 



(69) 



where r2jv(z, u) is the grand partition function for a ring 
of N sites. We then get 



(H d )/N = 2Prob(00000)\/i- 
+2Prob(10001)A 



4Prob(10000) % /5^ 
^ (70) 



Here Prob(OOOOO) is the probability of finding a randomly 
selected set of five consecutive sites all unoccupied in 
the equilibrium ensemble corresponding to the Hamilto- 
nian H' cl . These probabilities are also easily calculated. 
Treating z and u as variational parameters, we find that 
(Hd) is minimized for z — 0.35198 and u = 1.3752. For 
these values one finds that the density is p — 0.1952, 
Prob(OOOOO) = 0.28066, Prob(10000) = 0.082804, and 
Prob(10001) = 0.02443. These give 



E < -0.60333 



(71) 



which is an improvement over Eq. (|66[) , and quite close to 
the extrapolated value of —0.60356 obtained from Table 
I. This may be further improved by taking third-neighbor 
interactions in the classical Hamiltonian, but this will not 
attempted here. 



VI. STUDY OF GROUND STATES IN OTHER 
SECTORS 



We define a more general Hamiltonian 



H(\) = H Kit + \Y J W n 



(72) 



Since the W^s commute with Hku, all the eigenvectors 
of Hku can be chosen to be simultaneous eigenvectors 
of -ff(A), for all A. However, if we vary A, we can get 
different eigenvectors to have the lowest energy. 

Clearly, if A is large and positive, the ground state will 
lie in the sector with all W n = —1. Conversely, if A is 
large and negative, the ground state is the lowest energy 



A. Sectors with most W n 's positive 

Since the ground state for A = lies in the sector with 
all W n = +1, we have Af > 0. In fact, if the lowest 
excitation energy in the Hamiltonian Hkh is Ai?, we 
have Af = AE/2. At this point, the energy required to 
change a single W n from +1 to —1 becomes zero. We 
now study this sector using the variational techniques of 
Sec. V and try to estimate the difference between the 
ground state energy of this sector and the sector with all 

W n = l. 

Without loss of generality, we may assume that in this 
sector, Wn = —1, and the rest of the W's are +1. The 
basis vectors in this sector are of type \xU), or \Vy), 
where U and V are all possible strings of length N — 1 
obtainable from the string zzz . . . z of length N — 1, using 
the substitution rule zz — s- yx. Let \tp) be the eigenvector 
corresponding to the lowest eigenvalue of H in this sec- 
tor. It is easy to verify that (xU\H\xU') and (Vy\H\V'y) 
are negative, for all U and U' , and V and V. But, 
(Vy\H\xU) are positive. This implies that (xU\ip) and 
(xU'\tp) have the same sign for all U and U'. Similarly 
(Vy\ip) and {V'y\t/j) have the same sign for all V and V. 
This suggests a variational wave function of the form 



i 

71 



y/Prob(U)\xlI) - ^ y/PTob(V)\Vy) 



(73) 

Here Prob(t/) and Piob(V) are arbitrary functions, 
satisfying the constraint 



Y Prob(£/) = Y Prob(F) = 1. 



(74) 



Each configuration U is in one-to-one correspondence 
with the configurations of a nearest-neighbor-exclusion 
lattice gas on a linear chain of length (N — 2). Define a 
chain configuration as C = {rii}. We put 71$ = 1 if and 
only if there is a y in U in the position i + 1, otherwise 
rii = 0. Note that the last element of U cannot be a y. 
We specify C by a binary string of length (N — 2) . To go 
from C to U, we first add a single to the binary string 
of C at the right end, and then use the substitution rule 
10 — > yx. The remaining zeros in C are replaced by z's. 
Similarly, we specify V also by a binary string of length 
(N — 2), with x — >• 1, y, z —¥ 0, and as the leftmost 
element of the resulting string is always a zero, it may be 
deleted. 

As in the previous calculation, we construct a classi- 
cal Hamiltonian to variationally estimate the parameters 
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Prob(C). In this case, there is no translational symme- 
try, and in general, the lattice gas will have a non-trivial 
density profile. This is taken into account by making the 
activities of the lattice gas in the classical Hamiltonian 
site-dependent. We write 

N-3 N-2 

H% = +00 ^ n i n i+i - X! V^i- (75) 

i=l t=l 

The probability of each configuration C of the lattice gas 
is then given by 

Prob(C) = exp[-0H% (C)]/n N _ 2 ({z t }), (76) 

with Zi — exp(/3/ii), and QN-2({zi}) is the grand parti- 
tion function of the open chain of N — 2 sites. 

We note that the matrix H is unchanged under the 
space reflection mJV + l- i, and at the same time ex- 
changing x and y. This can be built into our eigenvector 
by assuming that if V are strings corresponding to lattice 
gas configurations C, we set 

Prob(F) = Prob(t/), (77) 

where U is the string corresponding the lattice gas con- 
figuration C T , the transpose of C. 

The rest of the calculation is done as before. By con- 
struction, we have 

(VVarlVVar) = 1. (78) 

It is straightforward to express (ip\H\ip) in terms of the 
marginal probabilities of the different local configurations 
of the lattice gas, remembering that there is no transla- 
tional invariance. For example, we get 

(iplS^Sfl^) = -Prob c (ni = 0). (79) 

In the simplest case, we work with only two parame- 
ters, and set Z\ = z' , and Zj = z for i 7^ 1. We would 
like to estimate the difference of the ground state en- 
ergy in this sector and the ground state over all sectors. 
These energies are of order N, and to cancel the lead- 
ing linear A^-dependence, we have to set z equal to the 
optimal value z* = 0.4045 to get the best energy value 
of Eq. We assume that N is large, so that only 

the term in the partition function corresponding to the 
largest eigenvalue is kept. Extremizing over z' we obtain 
z' = 0.2537, and for this value 

(VVar|#|^var) = -0.60057V + 0.1875. (80) 

This implies the following bound on the lowest eigenvalue 
in this sector 

E' Q < NE + A. (81) 

with A = 0.1875 providing an estimate of the energy 
gap between the ground state and the first excited state 
of the Hamiltonian. This estimate can be improved by 



adding more parameters in the variational wave function, 
or equivalently in the classical lattice gas Hamiltonian. 
A two parameter wave function would have z' and z" at 
the two opposite ends. Extremizing with respect to these 
parameters we find the energy gap to be 0.1642 with z' — 
0.2537 and z" = 0.6670. A four parameter wave function 
would have fugacities zi,Z2,zn-3 an d adjustable, 
and the rest of the Zj's set equal to z* . Table HTT1 shows 
the improvement in the value of the energy gap with the 
number of parameters used. 

We thus obtain a variational estimate of the energy 
gap of the first excited state from the ground state en- 
ergy. This matches quite well with the numerical esti- 
mates obtained in the previous section. 



Number of 


A 


parameters 




1 


0.18751 


2 


0.16419 


4 


0.15845 


6 


0.15642 


8 


0.15578 


10 


0.15556 



TABLE III: Estimate of the energy gap with the number of 
parameters used. 

It is straightforward to extend this treatment to sectors 
with two or more W n 's negative. There is an energy A 
required to create a single negative W n . Thus = A/2. 
If two defects are spaced far apart, the energy required 
to create two defects will be nearly 2A, with the cor- 
rection term decreasing exponentially with the distance 
between the defects. For n defects, the energy would be 
minimized if the defects are equally spaced. Thus the 
distance between the defects is N/n, and the energy cost 
of creating n defects AE(n) in H(X), for small n, is well 
approximated by 

AE(n) « - 2n\+ nA + nAexp(-BN/n), (82) 

where A and B are some constants. This then implies 
that for A = X\ + e, the density of defects in the true 
ground state of H(X) will vary as 1/| loge|. 

B. Sectors with most Wn's negative 

We now discuss the behavior of the ground state energy 
near the critical point X^,. This depends on the behavior 
of the ground state energy in sectors in which only a few 
of the WnS are +1. 

For a ring of N sites, the sector with all W n = — 1 
contains only two states, xxxxx ■ ■ ■ and yyyyyy ■ ■ ■ ■ The 
two are degenerate, with eigenvalue equal to —XN. 

Now consider the sector with only one W n = +1, say 
Wq = +1. Consider the state = \zxxxx ■ ■ ■) in this 
sector. From Eqs. (|5"3"|) , under Hkh, we have zx ^ yz, 
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and this state can make a transition only to the state 
"02 = \yzxxx ■ ■ ■) . And ip 2 can return to tp\ or go to 
-03 = \yyzxx ■ ■ ■). Thus, the dynamics may be considered 
as the dynamics of a particle z, which can hop to a nearest 
neighbor under the action of the Hamiltonian. There 
is a string of y's connecting the current position of the 
particle to the leftmost allowed position which is n = 1. 
This string can become longer, or shorter, as the particle 
moves, with no energy cost. When the z-spin is at the site 
N, it cannot move further to the right. The ground state 
energy Ef_ sector if this sector is seen to be the same as 
that of a particle with nearest-neighbor hopping, confined 
to move in the space 1 < x < N . It is thus given by 



R 9 

1 — sector 



-(N — 2)A - 2Jcos( 



N+l 



)• 



(83) 



Thus, we see that for large N, the state with all W n 's 
equal to +1 is no longer the ground state for A < J. 

We now consider a sector with exactly two of the W n 's 
equal to +1, and the rest negative. Let us start with 
the state \zxxx ■ ■ ■ zxxx • • • ), where the spins at two sites 
i = l and i = m + 1 < N are in the state z (these 
states will be referred to as z-spins in the following). This 
corresponds to Wn — W m = +1. Then, under the action 
of H(X), this state mixes with other states where the 
positions of the z-spins can change; the general state in 
this sector may be labeled by the positions of the z-spins, 
n and r 2 . We will write the vector as |t"i,?-2), where 
1 < n < m < r 2 < N . Then, for 1 < n < m and 
m + 1 < r 2 < N, we get 

H K it\ri,r 2 ) = -\n,r 2 + 1) - \n + l,r 2 ) 

-\r u r 2 -l}-\n -l.ra). (84) 

If the first z-spin is at m and the second is not at 
m + 1 , the first spin cannot move to m + 1 , as that site 
would be in spin state y, and the state zy cannot change 
[Eqs. ([52])], Similarly, if r 2 = N, and n ^ 1, then the 
second spin cannot move to the right. However, if the 
two z-spins are adjacent, then they can change to a state 
zz =± yx [Eqs. (|52l) ]. But from the state yx the state can 
only return to zz. 

If we disallow the transitions to state yx, the z-spins 
act as independent particles moving in two disjoint re- 
gions of space, 1 < r\ < m and m + 1 < r 2 < N. In 
this case, the minimum energy of this system is just the 
sum of the energies of two particles. This energy is an 
upper bound on the true ground state energy of this sys- 
tem. Thus, we find that the ground state energy in this 
sector, E2_ sector , has the upper bound 



ES 



2 — sector 



< -2Jcos( 



m + 1 ' 



2Jcos( 



N-m + 1 



-X(N-4). 



(85) 



Next, suppose that the state with the m-th site in the 
y-state and the (m + l)-th in the cc-state is called the 
state n = m + 1, r 2 = 1, and a similar definition for the 
other end. Then the range of r± is at most m + 1, and 



the range of r 2 is at most I — m + 1. By excluding some 
states (here r\ = m + 1, r% ^ r*i), the kinetic energy can 
only increase, and hence we have 

ELsector > -2JC0s(-^-)-2JcOs(— ^— — ) 



m 

-A(iV-4). 



' N - m + 2 ' 



(86) 

For JV,m > 1, these bounds can be expanded in pow- 
ers of 1/m, and have the same leading order correction. 
Also, the minimum energy corresponds to equally spaced 
defects, with m = N/2. 

We can easily extend the discussion to sectors with 
three, four or more W n 's equal to +1. In case the 
lengths of the intervals between the positive W n 's are 
TOi, 7T7 2 , 777,3, ' ' * , m r , the bounds on the lowest energy in 
this sector E a r _ sector become 



2J£, 



777i + 2 



) - X(N - 2r) < E 9 r _ 



V 

<-2jVcos( — - — )-X(N-2r). (87) 

^ TOj + 1 

1=1 

Thus, we see that for A > J, the ground state belongs 
to the sector with all W n '8 equal to — 1. If A = J(l — e), 
the ground state will be in the sector with 77 equispaced 
bonds with W n — +1, where the spacing I between them 
w A^/77 is given by e -1 ^ 2 . The minimum energy per site 
of H (A) for A = J(l — e) varies as e 3 / 2 for small e. Equiv- 
alently, if we restrict ourselves to sectors with only a frac- 
tion e of W n 's having the value +1, the minimum energy 
per site varies as — e 3 / 2 . This is equivalent to the state- 
ment that for Hku corresponding to A = 0, in the sector 
with the fractional number of positive Wn's being equal 
to A, the minimum energy per site varies as A 3 / 2 . 



VII. DISCUSSION 

In this paper, we first analyzed the symmetries of a 
spin-5 Kitaev chain. We found a Z 2 invariant, W n , asso- 
ciated with every link (77, 77 + 1), namely, invariants for 
the model defined on a ring with N sites. For integer S, 
these invariants commute with each other and the Hamil- 
tonian. The Hilbert space can therefore be split into 2 N 
sectors, where the Hamiltonian is block diagonal. For 
half-odd integer S, W n anti-commutes with W n ±i and 
commutes with the rest. We showed that this implies 
that all the eigenstates of the half-odd-integer spin mod- 
els are 2 N / 2 -io\d degenerate, thus showing a qualitative 
difference between the integer and half-odd-integer mod- 
els. We have developed a formalism to compute the di- 
mensions of the invariant sectors. We showed that the di- 
mension of most of the sectors can be calculated in terms 
of products of 2 x 2 matrices T + and T_ . For 5=1 the 
quantum dimension of the sector with all W n — 1 is the 
golden ratio, (1 + \/5)/2. For S — > 00, the quantum di- 
mension tends to S + 1/2 in both the W n = 1 and the 
W n = — 1 sectors. 
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We have then studied the spin-1 case in detail. We have 
found that the ground state lies in a sector which can be 
mapped to a quantum lattice gas model with nearest- 
neighbor exclusion. We developed a variational wave 
function that relates the quantum mechanical averages 
to the correlation functions of a classical lattice gas with 
nearest-neighbor exclusion. We considered a more gen- 
eral Hamiltonian with a term proportional to the sum of 
the conserved quantities, and showed that as a function 
of the coupling constant A, this would show gapless exci- 
tations in the range A^ < A < A| . We extended our vari- 
ational calculation to study how the ground state energy 
and the defect density would vary near the two critical 
points AJ and A§. At A = Af , Eq. (|82|) implies that the 
energy of the lowest excited state in a system of length 
L goes as E ~ exp(— BL), corresponding to a state in 
which one W n = — 1 while all the other W n = 1. By the 
usual scaling arguments, the gap to the first excited state 
goes as 1/ L z , where z is the dynamical critical exponent. 
We therefore conclude that z — oo. At A = A§, the low- 
energy excitations form a low-density gas of hard-core 
particles. In one dimension, this can be mapped to a 
system of non-interacting spinless fermions with a non- 
relativistic spectrum E ~ k 2 . Hence in a system of size 
L, the gap to the lowest energy states goes as 1/L 2 corre- 
sponding to k ~ thus z — 2. It would be interesting 
to find the value of z in the critical region A J < A < 

Finally, we note that there is another interesting one- 
dimensional spin model called the golden or Fibonacci 
chain [38l . [39J , for which the number of states on a ring 
of size N is the same as that of the spin-1 Kitaev chain 
in the sector with all W n = 1. The Hamiltonian for this 
model is 

H GC = + "i-i - 1) 

i 

-n l _m l+1 ( 7 - 3 / 2 ( Tf + 7" 3 n t + 1 + 7 2 )) , 

(88) 

where rii = (1 — er z )/2. It has been shown [H, that 
this model is critical. Its long-range correlations are de- 



scribed by a SU(2) level 3 Wess-Zumino-Witten model, 
which is a conformally invariant field theory with cen- 
tral charge equal to 7/10. The Hamiltonian in Eq. (|88|) 
differs from the spin-1 Kitaev chain in the W n — 1 sec- 
tor by terms which are products of the n% operators. We 
have shown that the spin-1 Kitaev chain is gapped. Thus 
these terms correspond to some relevant operators which 
take the golden chain Hamiltonian away from criticality. 

We can show that it is possible to add multi-spin terms 
to the minimal Kitaev chain which reduce to the extra 
terms in the W n = 1 sector. We need to add products of 
the rii operators to the minimal Kitaev chain to obtain 
the golden chain in the sector with all W n = 1. The basis 
states |t), | I) that we use in Sec. [IV] are eigenstates of 
the rii operators with eigenvalues 1 and respectively. 
The | t) state represents a state with the head, namely 
\x). The | \) state represents either an empty site, \z), 
or a tail, \y). It is clear from Eqs. (|47]l and (|l8| that 
the operator P x = (1 + Yj x )/2 has eigenvalues 1 for \x) 
and for \y) and \z). Since all the £ matrices commute 
for integer S, they commute with the invariants and are 
block diagonal within the invariant sectors. Thus, the 
Hamiltonian, 

Hkgc = 7~ 3/2 ^c - £ (l - P? +1 - 

i 

+7 a i?_ii?+i + - :i i'f :i';r;. : ) , (89) 

when restricted to the W n = 1 sector, is exactly the 
golden chain Hamiltonian discussed by Feiguin et at and 
others [H, [3§| . We have thus constructed a realization 
of the golden chain model as a spin-1 chain. 
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